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Abstract. We investigate a state estimation problem for the dynamical 
system described by uncertain linear operator equation in Hilbert space. The 
uncertainty is supposed to admit a set-membership description. We present 
explicit expressions for linear minimax estimation and error provided that 
any pair of uncertain parameters belongs to the quadratic bounding set. We 
introduce a new notion of minimax directional observability and index of non- 
causality for linear noncausal DAEs. Application of these notions to the state 
estimation problem for linear uncertain noncausal DAEs allows to derive 
new minimax recursive estimator for both continuous and discrete time. We 
illustrate the benefits of non-causality of the plant applying our approach 
to scalar nonlinear set-membership state estimation problem. Numerical 
example is presented. 

Key words, set-membership state estimation, minimax, uncertain linear 
equation, DAE, descriptor systems, implicit systems, Kalman filter. 

1 Introduction and problem statement 

The applications of differential-algebraic equations (DAEs or descriptor sys- 
tems) in economics, demography, mechanics and engineering are well known [1] 
This in turns motivates researchers to investigate DAEs from the mathemat- 
ical point of view [2]. Here we focus on a design of state estimation algorithm 
for uncertain linear non-causal DAE. 

The most common approach to DAEs investigation is to reduce it to some 
canonical form which in turn is equal to some normal ODE. In particular, 
one of the basic results of the algebraic theory of regular linear DAEs with 
constant matrices 1 was introduced in [3]: if the linear DAE with constant 

^erc wc present a formulation for the finite-dimensional Banach space. 
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matrices 



Fx = Cx + Bf 



(1) 



is well defined (det[AF — C] ^ ) then for all initial values a; (to) = there 
exists the unique solution x(-) provided that /(•) is sufficiently smooth. The 
index s of the pencil F,C is said to be an index of linear DAE (1). One 
can reduce (1) to the ODE via change of coordinates so that the pencil F, C 
brings into canonical form [1] and differentiating exactly s times provided 
that / is sufficiently smooth. In such a way one can derive an analogue of 
the celebrated Cauchy formula for the linear regular DAEs with constant 
matrices. This result is generalized to variable coefficients by means of a 
standard canonical form (SCF): in [5] it was shown that (1) with analytical 
F, C, B is solvable (i.e. for every sufficiently smooth / there exists at least 
one continuously differentiable solution to (1) provided F, C to be sufficiently 
smooth) if there exists the SCF for (1). Note that in this case rankF(t) 
changes only at finite number of points from within any compact [to?^1- 
In [6] it was noted that not all solvable DAEs can be put into SCF and the 
solvable DAE is equal to some differential-algebraic equation in the canonical 
form which generalize SCF. In this respect we say that DAE is causal if it 
can be reduced - at least locally in nonlinear case - into normal ODE. The 
geometry of the reduction procedure for nonlinear causal DAEs F(x, x) — 
was investigated in [7, 8] , where the index of DAE was defined as a smallest 
natural s so that the sequence of the constraint manifolds [7] 



becomes stationary for k > s. This coincides with the definition of the index 
of linear DAE. Further discussion of the DAEs solvability theory and related 
topics is presented in [1, 2]. 

The noncausal DAE differs radically from the causal one. For instance, 
consider 



Let a^C') G I^to,^), /(•) 6 Li(io,T) and xq G K. By definition put 



M k := TW k -i D M fc _!, M := {(x,p) : F(x,p) = 0} 
W := {x G R n : (x,p) G M } 



"](£)= [ Cl C2 ](S)+/(*),b o](5)(* ) = a:o 



(2) 



xi(t) := exp(ci(t - t ))x + / exp(ci(t - s))c 2 x 2 (s) + f(s)ds 
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It is clear that any solution of (2) is given by the formula 1 1— > (xi(t), X2(t)) T . 
According to a behavioral approach [9] one can think about x 2 as an input 
or as a part of the system state representing uncertain inner disturbance 
generated by the plant itself. In order to clarify this ambiguity we shall give 
an exact definition of the DAEs solution accepted in this paper. According 
to [10] x(-) is said to be a solution of 

jFx{t) = C(t)x(t) + f(t) (3) 

with initial condition Fx(to) = if Fx(-) is totally continuous function, x(-) 
satisfies (3) almost everywhere and Fx(to) = holds. This definition allows 
to properly define the adjoint system. Also it guarantees that the linear 
mapping induced by (3) is closed in corresponding Hilbert space [10]. Note 
that this doesn't hold for (1). Another useful application of the introduced 
solution is in the control theory. In [11] authors discuss difficulties arising 
while applying of proportional feedback / = Kx to the (1): even well defined 
DAE (det(sF-C) ^ 0) may become singular (det(sF-C-BK) = 0). In [12] 

a properly stated leading term A(t)—F(t)x is used in order to give a feedback 

solution to LQ-control problem with DAE constraints. This generalizes the 
definition of DAE solution [10] to the case of variable matrices. 

Recently solvability conditions for abstract semi-linear non-causal DAE 
has been studied in [13] assuming that the pencil sF — C is singular, F, C 
are closed linear mappings in abstract Banach space. Properties of the solu- 
tions of noncausal implicit differential equation with special structure were 
discussed in [14]. Note that non-causal DAEs are not just a "pure" mathe- 
matical structure which is suitable for solving control or observation problems 
only - some potential applications of non-causal DAEs was briefly discussed 
in [15]. 

A state estimation framework for linear dynamic models has several widely- 
used approaches: H 2 / filtering and set-membership state estimation. H 2 - 
estimators like Kalman or Wiener filters (also known as minimum variance 
filters [16]) give estimations of the system state with minimum error vari- 
ance. These filters require an exact model of signal generating process and 
full information about a statistical nature of noise sources. Recently, the 
^-estimation for linear DAEs has been studied in [17]. Authors derive a 
so-called "3-block" form for the optimal filter and a corresponding 3-block 
Riccati equation using the maximum likelihood approach. The obtained re- 
cursion is stated in terms of a block matrix pseudoinverse. In [18] the filter 
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recursion is represented in terms of a deterministic data fitting problem solu- 
tion. Authors introduce an explicit form of the 3-block matrix pseudoinverse 
for a descriptor model with special structure, so that the form of obtained in 
[18] filter coincides with presented in [17]. A brief overview of steady-state 
if 2 -estimators is presented in [19]. Optimal estimators minimize the 2- 
induced norm 2 of the operator that maps unknown disturbances with finite 
energy to filtered errors [20]. In literature it is common to construct subop- 
timal estimators [21] that guarantee aforementioned norm to be less then a 
prescribed performance level 7. Note that estimators are certain Krein 
space H 2 filters [22]. Krein space approach was used in [23] for risk-sensitive 
filtering in linear time-invariant (LTI) descriptor models with regular matrix 
pencil under stochastic noise. A linear matrix inequality approach was used 
in [24] in order to construct reduced order ifoo-filter for LTI DAE with reg- 
ular matrix pencil. An up to date description of the state of the art is to be 
found at [2 5]. 

In the sequel we focus on the following problem: given some element 
(for instance measurements of the system output) y from some functional 
space one needs to estimate the expression £{9) provided that g(9) = 0. This 
problem becomes non-trivial if the latter equation has more than one solution 
and the equality y = C{9) holds. In this case the estimation problem may be 
reformulated as follows: given y = C(9),9£<d,y£Y one needs to find the 
estimation £{9) of the expression £{9) provided that g{9) = and C (•),£(•) 

are given functions. Note that £{9) := £{9) if the equation y = C(9) has the 
unique solution 9. 

The estimation problem is said to be linear if 0, Y are linear spaces and 
C (■),£(■) are linear mappings. It is a common case when 

C{9) = H<p + Dr,, g{9) = L<p + Bf, (*) 

where 9 = (x,f,rj) C X x F x Y, H,D,L,B are linear mappings. The 
linear estimation problem is said to be uncertain if D ^ 0, L and B are 
non-trivial or if B = and N(L) = {if : Lip = 0} ^ {0}. Note that the 
choice of the solution method depends on the "type of uncertainty" : if /, 77 
denote realizations of random elements then it's natural to apply probability 
methods. This requires an a priori knowledge of distribution characteristics 
of the random elements. In the sequel we assume that there is uncertainty 

2 sup of the relation between the Euclidean norms of estimation errors and model dis- 
turbances 
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in (*) if distributions of random elements or some deterministic parameters 
of the system are partially unknown. It is natural to choose the estimation 
from some class to be optimal in the sense of the given criteria. According 
to this the linear uncertain estimation problem is said to be minimax if 
the class of estimations £(ip) is restricted to all linear functions (u, y) + c 
of y and the criteria is set to be the minimum of the worst-case error. A 
description of the state of the art in the theory of linear uncertain minimax 
estimation problems with special £, L, H, B, D in special spaces is to be found 
at [26, 27, 28, 29, 30, 31, 32, 33, 34]. 

1.0.1 Author's own research activities 

Classical theory of uncertain estimation problems [26]- [29] works well when 
the linear mapping L in (*) has bounded inverse. One of the author's theo- 
retical achievements is the extension of the linear minimax estimation theory 
to abstract equations with closed linear non-injective mapping [35, 36, 37, 38] 
in Hilbert space. This extension is based on the general duality principle as- 
serting that the linear minimax estimation problem is equal to some control 
problem with convex non-smooth cost and linear constraints provided that 
uncertain parameters belong to closed bounded convex sets in corresponding 
Hilbert spaces [39, 40, 36]. Note, that these results were previously obtained 
for the finite dimensional space [41, 42, 43, 44, 44, 45]. 

In order to apply the abstract theory to DAEs the sufficient conditions on 
DAEs matrices were introduced asserting that the linear mapping induced 
by the noncausal DAE is closed and has closed range [10]. Also a general- 
ization of the integration by parts formula and the necessary and sufficient 
conditions of solvability of DAE in the form (3) is presented in [10]. The 
solvability condition is obtained via application of Tikhonov regularization 
approach. With help of these results new notions of the minimax directional 
observability and index of causality for discrete time [16] and continuous 
time linear non-causal DAEs [47] were introduced. The minimax directional 
observability provides a qualitative description of DAEs singularity with a 
respect to the given observations. Using this notion author developed sev- 
eral representations of the minimax estimation [48, 49, 50]. The final result is 
an algorithm which allows to compute the minimax estimation of the linear 
noncausal DAE state in the real time [47]. The structure of this algorithm 
coincides with celebrated Kalman filter recursions for normal linear ODEs 
with continuous time. Similar results were obtained for linear noncausal 
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DAEs with discrete time [51, 52, 53, 54, 55, 56, 46]. In [46] the author gives 
a complete solution to the problem of recursive implementation of the mini- 
max a-posteriori estimation (similar to posed in [28]) for the linear non-causal 
DAEs with discrete time. 

Also the theory developed in [36] was applied to the state estimation of 
the solutions of finite-dimensional linear boundary-value problems [57] and 
the minimax mean-square estimations of trends [58]. 

Notation. c(G, •) = sup{(z, /), / G G}, 
S(G, x) = if x G G and +oo otherwise, 
dom/ = {x G H : f(x) < oo}, 

f*(x*) = sup x {(x*,x) - f(x)}, (L*c)(u) = M{c(G,z),L*z = u}, 
(fL)(x) = f(Lx), (L*c)(u) = mf{c(G,z),L*z = u}, 
cLf = /**, ^ 

Argini u f(u) denotes the set of minimum points of /, 
df(x) denotes the sub-differential of / at x, 
(•, •) denotes the inner product in Hilbert space, 
S > means (Sx, x) > for all x, 
L* denotes adjoint operator, 

Pl* denotes the orthogonal projector onto R(L*), 

R(L), N(L) and @(L) denote the range, the null-space and the domain of 

the linear mapping L, 

F' denotes transposed matrix, 

F + denotes pseudoinverse matrix, 

E denotes the identity matrix, 

diag(Ax . . . A n ) denotes diagonal matrix with A i; i — 1, n on its diagonal, 
G denotes the closure of the set G, 

[x\, . . . ,x n ] denotes an element of the Cartesian product H 1 x ■ ■ • x H n of 
Hilbert spaces Hi,i = l,n, 

MJ 1 denotes n-dimensional arithmetic Hilbert space, 

C mxn (t , T) denotes the space of all continuous on (t , T) functions with val- 
ues in R mxn , 

h2(to,T) denotes the space of all measurable functions with finite integral 

W'z'ftojT) denotes the space of all absolutely continuous functions with 
derivative from L^tcb T), 

M£ denotes expected value of the random vector £. 
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2 Linear uncertain estimation problem 



In this section we present the main result of the paper [36]. All proofs are 
given in [36]. 

Suppose that L(p G 3f and 

y = H V + r] (4) 

The mappings L, H and the set £f are supposed to be given. The element 
r\ is uncertain. Our aim is to solve the inverse problem: to construct the 
operator mapping the given y into the estimation £(ip) of expression £(ip) and 
to calculate the estimation error a. Now let us introduce some definitions. 

The operator L : H \— > T is assumed to be closed. Its domain @(L) is 
supposed to be a dense subset of the Hilbert space H, H G Sf(7{,y). Note 
that the condition Lip G £f is equal to the following 

L<p = f, (5) 

where / is uncertain and belongs to the given subset of the Hilbert space 
T. In the sequel 77 is supposed to be a random 3^-vahied vector with zero 
mean so that its correlation R n G where 3% is some subset of J?(y,y). 
Also we deal with deterministic 77 so that (/, 77) G £7, where Q is some subset of 
J 7 xy. Note that the realization of y depends on 77, H and /. Also it depends 
on elements of N(L) — {<p G ^(-^) : -^V 9 = 0} so that y = H(ipo + tp) + 77, 
where y?o may be thought as inner noise in the state model (5). 

Let £(<f) = (£, cp), £(cp) = (u, y) + c. Since L, H are not supposed to have 
a bounded inverse mappings the £(<p) and £{ip) are not stable with a respect 
to small deviations in /, 77. Also /, 77 are supposed to be uncertain. Therefore 
we use the minimax design in order to construct the estimation. 

Definition 1. The function £(<p) = (u, •) + c is called the a priori minimax 
mean-squared estimation iff o~(£,u) = inf U)C a(£, u) where 

a(£,u):= sup M(%)-%)) 2 (6) 

The number a(£) = a^(£,u) is said to be the minimax mean-squared error 
in the direction £. 
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On the other hand the a posteriori estimation describes the evolution of 
the central point of the system reachability set 

(Lip, y - Hp) G G 

consistent with measured output y [26, 28, 27]. Note that the condition 
(L(p,y — Hep) G Q holds if \\y\\ < C for some real C. But it doesn't hold in 
our assumptions if t] is random since ||-R r /|| < c doesn't imply |||/|| < C for 
realizations of 77. Therefore 77 is supposed to be deterministic. 

Definition 2. The set 

Xy = {<pe 9(L) : (Lip, y-Hp)e Q} 

is called an a posteriori set. The vector (p is said to be minimax a posteriori 
estimation of tp in the direction i (f-minimax estimation) iff 

d(£) := inf sup \(£,<p) - "0)1 = sup \(£,<p) - (£,ip)\ 

The expression d(t) is called the minimax a posteriori error in the direction 
I [I- minimax error). 

In the sequel the minimax mean-squared a priori estimation (error) is 
referred as minimax estimation (error). 

Proposition 1. Assume that , & are convex bounded closed subsets of J 7 , 
■^{y 1 y) respectively. For the given I e 7i the minimax error &(£) is finite 
iff 

£ - H*u e domcl(L*c) n (-l)domcl(L*c) (7) 
for some u G y . Under this condition 

a(£,u) = sup (R v u,u) + 

(8) 

- [c\(L*c) {£ - H*u) + cl(L*c) {-£ + H *u)f 

where 

R(L*) C domcl(L*c) C R(L*) 
If Arginf u cr(£, u) ^ 0, then £(p) = {u, y) + c, where 

u G Arginf u <r(£, u) 
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and 

c = ^(cl(L*c)(£ - H*u) - c\(L*c)(-£ + H*u)) 

Theorem 1. Suppose that <S is convex bounded closed balanced set and G 
int£f. Also assume that 

7] G {q : M(rj,rj) < 1} 

Then for the given I G 7i the minimax estimation a(£) is finite iff I — H*u G 
R(L*) for some u G y '. Under this condition there exists a unique minimax 
estimation u and 



o~(£, u) = min a(£, u) 

u 

a(£, u) = (u, u) + min{c 2 (^, z), L* z = £ - H*u} 



(9) 



If R(L) , H (N (L)) are closed sets then u is determined by the following con- 
ditions 

u - H Po G H(dI 2 (H*u)), Lp = 0, 

I 2 {w) = mm{c 2 {&,z),L*z = P L *(£-w)}, ^ 

z 

Corollary 1. Let 

&={feF:(fJ)<l}, V e{i 1 : Mfari) < 1}, 
and suppose that 

1) R(L),H(N(L)) are closed sets; 

2) R(T) = {[Lx, Hx],x G @{L)} is closed set. 

Then the unique minimax estimation u is given by u = Hp provided that 
£ G R(L*) + R(H*), p obeys 

L*z = £- H*Hp, , 

11 

Lp = z V 1 

The minimax error is given by the following expression 
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Corollary 2. Assume that linear mappings L : TC > T, H G Jz?(7Y, y) 
obey 1) or 2) (Cor. 1). Then (11) has a solution z G @(L*),p G @{L) iff 
£ = L*z + H*u for some z G @(L*), u ey. 

Corollary 3. Under conditions of Cor. 1 for any £ G R(L*) + R(H*) and 
some realization of y(-) we have (u,y) = {£,0), where (p obeys 

L*q = H*(y-H0), 
L<p = q 

Consider an a posteriori estimation. 
Proposition 2. Let Q be a convex closed bounded subset ofyxj 7 . Then 
R(L*) + R(H*) C domc{X y , •) n (-l)domc(^, ■) C 

(13) 



R(L*) + R(H*) 

The minimax a posteriori error in the direction £ is finite iff ' £ G domc(X y , -)n 
(— l)domc(X y , •) and 

(l,0) = hc(Xy,£)-c(Xy,-£)), 

i (14) 
d(£) = -(c(X y ,£) + c(X U7 -£)) 

Theorem 2. Let 

e = {(/,77):||/H 2 +N| 2 <l}, 
and assume that 1) or 2) from Corollary 1 holds. The minimax a posteriori 
estimation (p obeys 

L*q = H*(y-H0), 
Lip = q 

iff £ G R(L*) + R(H*). The estimation error is given by 

d(£) = (l-(y,y-H0))h(£) (16) 

Corollary 4. Assume that the conditions of Theorem 2 are fulfilled and 
£(<p) = (£,0) for any £, where (p obeys (15). Then 

inf sup \\<p — x\\ = 

sup \\<p — x\\ — (1 — (y, y — H0)a max a(£) 

XdXy ll^ll = i 
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In order to apply these results for linear DAEs we investigate some prop- 
erties of the linear mapping induced by DAE [10]. Let 

bp(t) = [±F<p(t)-C(tMt),F<pfo)] 

and set 

®(L) = W F := M-) G L 2 (t , T):t^ F<p(t) G W 2 (t , T)} 
It is clear that L<p(t) = [f(t), fo] is equal to 

jFip{t) = C{t)<p{t) + f{t),F<p{t ) = fo 

Next proposition describes the adjoint L*. 
Theorem 3. If x(-) G W F , z G W F > then 

j\j t Fx(t),z(t)) + (jFz(t),x(t))dt = 

(Fx(T),F' + Fz(T)) - (Fx(t ),F'+F'z(t )) 

L is closed linear mapping and its adjoint L* : h 2 (to,T) xK"-> L 2 (to 5 ^) is 
defined as follows 

L%z,z )(t) = -j t F'z(t)-C'(t)z(t), 

@(L*) = {(z, F l+ F'z{t Q ) + d): ze W F ; F'z{T) = 0, F'd = 0} 

Note that R(L) is not necessary close. A sufficient condition for R(L) to 
be close is introduced in the next theorem assuming 3 that 

f=(« »).c=(a g) 

Theorem 4. If 

sup \\Q(e)C' 2 \\ mod < +oo, Q(e) : = (e 2 E + C' 4 C,)-\ \\F\\ mod := ^ \F t 



(17) 



1>E>-1 

hJ 



then R(L) is closed. 

Next subsections demonstrate the application of the above theory to the 
linear estimation problem for linear DAEs. 

3 This assumption holds for any linear DAE with constant matrices. 
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2.1 DAEs with continuous time 



In this subsection we present the main result of preprint [47] - linear reduced 
order minimax filter for linear noncausal DAEs with continuous time. All 
proofs are given in [47]. 

Consider a pair of systems 

jFxit) = C(t)x(t) + f(t), Fx(t Q ) = 0, 

y(t) = H{t)x(t)+ri{t),te[t ,T\, 

where x{t) G W 1 , f(t) G K m , y(t) G W, rj(t) G W represent the state, 
input, measurement output and measurement noise respectively, F G M mxn , 
t i * C(t) G C mx "(t ,T), /(•) G L 2 (t , T),t ^ H{t) G CP x "(t ,T), t ,T G R. 

According to [10] we say that x(-) is a solution of (18) if Fx(-) G W™(to,T), 
the derivative of Fx(-) coincides with the right side of (18) almost everywhere 
and Fx(to) = holds. 

In the sequel we assume that rj(-) is a realization of the random process 
rj with zero mean satisfying 

rj e W = {77 : M f (R(t)r,(t),r,(t)) < 1} (19) 

J to 

and 

/(•) G G = {/(-) : f T (Q(t)f(t)J(t))<l}, 

J t 

where Q(t) G M m x m , Q = Q' > 0, i?(t) G M pxp , R' = R > and Q(t),R{t), 
i? _1 (t), Q _1 (t) are continuous functions of i on [to,^ 1 ]- 

Suppose y(t) is observed in (18) for some x(-), f G G and 77. Our aim here 
is to construct an algorithm 4 giving online estimation of the linear function 

x(-) ^ (£,Fx(T)) 

on the basis of the measured on [to,^ 1 ] realization of the output y{t). With 
this purpose we introduce a notion of the linear minimax estimation [36]. 

Definition 1. The function u(y) = J^(u(t), y(t))dt is called minimax mean- 
squared a priori estimation if 

inf a{u) = a(u) 

u 

4 In literature it is common to refer to this algorithm as filter [59] 
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where a{u) = s\ip x /\j(.\ v M[(£,Fx(T)) — u(y)} 2 is a maximum estimation 
error for «(•). The number a = cr(tt) is called a minimax mean-squared a 
priori error. The state x(t) is said to be minimax observable in the direction 
£ iff a < +00. 

The minimax directional observability differs from the classical observ- 
ability property in the following way. If system state x(s) is minimax observ- 
able in the direction I then the projection of the reachability set (consistent 
with some realization of y(t),t < t < s) onto direction I is expected to 
be [— a, a], where a denotes the minimax estimation error. The expected 
estimation error varies in [0, a] and depends on the noise realization, ini- 
tial condition and input. If x(s) is unobservable in the minimax sense for £ 
then the minimax estimation is set to zero and b = +00. This means that 
the structure of the measurements do not provide any information about 
(£, x(s)). Therefore the minimax directional observability provides a qualita- 
tive description of DAEs singularity with a respect to the given observations. 
In particular, regular DAE is observable in the minimax sense for any direc- 
tion in contrast to the classical observability. 

Remark 1. This definition generalizes the notion of linear minimax a priori 
estimation introduced in [27]. Here we follow a common way of [27] deriving 
the minimax estimation: first step is to describe a dual control problem, next 
step is to solve it and the last step is to derive a minimax filter. 

Assume that u(-) E h 2 (t ,T), £ G M m and z(-) obeys DAE 

jF'z{t) = -C'{t)z{t) + H'(t)u(t), F'z(T) = F'l (20) 

Let v(-) denotes any solution of homogeneous DAE (20). Next proposition 
gives a generalization of the celebrated Kalman duality principle [59] . 

Proposition 1. The minimax estimation error 

a{u) = sup M[(£, Fx{T)) - u(y)} 2 -> inf 

is finite iff (20) has a solution z(-). The minimax estimation problem o~(u) — > 
ini u is equal to the following optimal control problem 

I(u) = min{/ (Q~ 1 (z-v),z-v)dt}+ f (R^u, u)dt -»• min, (21) 

v J to J to u 

provided that z(-) is some solution of (20). 
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Proposition 1 states that minimax estimation problem is equal to some 
optimal control problem for appropriate I which is called dual control prob- 
lem. In the next proposition we introduce a representation for the minimax 
estimation and error. 

Proposition 2. Let p(-) denotes some solution of the two-point boundary 
value problem 

^-Fx{t) = C(t)x(t) + Q-\t)z(t), Fx(t ) = 0, 
j (22) 
^-F'z{t) = -C'{t)z(t) + H'(t)R{t)H{t)p{t),F'z{T) = F'i 

Then minimax estimation u is given by u = RHp, the minimax error is 
represented as a = (£,Fp(T)). 

It is known that (18) may be converted into SVD coordinate system [60] 
so that without loss of generality we assume that 

f = ( e o 8),c=(g c d)Mt) = ( Q Q l = (g = ID 

where S = H'RH. By definition, put A(t) = C x - Q 2 Ql 1 C z - (C 2 - 
Q2g 4 " 1 C , 4)^(S , 3 + Qg 4 - 1 C , 3), M(t) = S 1 + C' 3 Q^C 3 -(S 2 C' 3 Q^ 1 C 4 )Si(S 3 + 
C' A Q:y 3 ), C = S+((Cz - Q 2 Ql 1 C±) T - (S 3 + C' 4 Q?C 3 )K), G(t) = Qi - 

Q2Q1 Qi + (C*2 — Q2Q4 Ci)S2{C2 — Q2Q4 c*4) T - 

Theorem 5. Assume that t h-> §£(t) = (S 4 (t) + C'^ifyQ^ 1 (t)C 4 (t)) is mea- 
surable matrix-valued function. For any f 6 I" the minimax estimation of 
the inner product (£,Fx(T)) is given by 

(£^x~{f)) = (£ 1 ,x(T)) 
where x is the solution of the initial-value problem 

jx = (A(t) - K{t)M{t))x + K(t)[E,C']H'Ry(t),x(t ) = 0, 
K = AK + KA' + KMK - G, K(t ) = 

The minimax estimation error is given by a = (£1, K{T)l\), where t is splitted 
into according to the block structure of F. 
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2.2 DAEs with discrete time 



In this subsection we present the main result of the preprint [16] - linear 
recursive minimax filter for linear noncausal DAE with discrete time. All 
proofs are given in [46]. 
Consider the model 

F k+1 x k+1 - C k x k = f k , F x = q, (24) 
y k = H k x k + g k ,k = 0, 1, . . . (25) 

where x k G M", f k G M m , y k , g k G MP represent the state, input, measurement 
output and measurement noise respectively, F k ,C k G M mxn , H k G MP xn and 
initial state Xq belongs to the affine set {x : Fqx = q}, q G M. m . In what 
follows we assume that 

£ e = <l,VneN} (26) 

where £ = [q, {f s }, {g s }}, {f s } = [f , f 2 , . . .}, £ k = [q, {f s }%, {g s }%}, {/ s }g = 
[fo ■ ■ ■ fk] is the projection of {f s } onto linear span of e\ . . . e k+ i, e\ = [1,0...], 
*n(^) = (Sq,q) + J2o~\S t f t J t ) + (Rigugi), S, S k G M mxm , R k G W xp are 
symmetric and positive-definite. 

Suppose that y\ is being observed in (25) with x k = x* k and g k = g* k 
provided that x* k is derived from (24) with f k = f£, q = q* and ^* = 
[q*j {fs}' idt}] e y • Our aim here is to describe the evolution in r of the 
^-minimax estimation x T of the state x* at instant k = r along with error 
p(£, r) through dynamic recurrence-type relation and to efficiently describe 
the structure of the minimax observable subspace C(r). For this purpose we 
shall apply the theory developed in the previous section. 

Definition 2. The set C(t) = {£ : p(£,r) < oo} is called a minimax ob- 
servable subspace for the model (24) at the instant k = r. Its co-dimension 
I T = n — rankQ r is called an index of non- causality of the model (24). 

Theorem 6. The minimax observable subspace for the model (24) is given 
by £(r) = {i : P+PJ = £} and 

X(t) = {xe R n : (P T (x -x T ),x- x T ) < /3 r } (27) 

where x T = P^r T , j3 T := 1 — a T + (P T x T , x T ) 

Pk = H' k R k H k + F k [S k ~i — Sk-iCk-iB^^C^^Sk-^Fk, 
Pq — FqSF + HqR H , B k = P k + C' k S k C k , 
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ati = aii-i + (RiHuyi) - (B^r^, rvi), a = (R Q y Q ,y ) and 

r k = i^5 , fc _iC fc _i J B+_ 1 r fc _i + H' k R k y k ,r = H'^R^yo 

If £ £ >C(r) then x T is l-minimax estimation of x* and p(£, r) = j3r {Pr~£, £)^ ■ 

Example. Consider the following filtration problem: given measure- 
ments y k one needs to construct the estimation x T of the state x k at instant 
k = t and to describe the estimation error provided that 

x k +i = c k x k + v k (x k ) + f k ,x = q,y k = h k x k + w k (28) 

V l— * v k{p) is some real- valued function and uncertain scalar parameters q, 
w k and f k are restricted by the inequality 

T-l 

Sq 2 + RkW 2 k + S k f k <l,S,R k ,S k >0 (29) 

k=0 

Let us show how one can construct x T by means of the set-membership state 
estimation approach for linear non-causal descriptor systems described in 
this section. Let z k = [zx tk ,Z2 t k\ obeys 

Fz k+1 = C k z k + f k , Fz = q, a k = H k z k + w k (30) 

where F = (1, 0), C k = (c k , 1), H k = (h k , 0). Note that for any real z 2 , k there 
exists exactly one z± jk so that z k obeys the first equation in (30). 
We shall apply Theorem 6 in order to construct the £-minimax estimation 
of z T . Using definitions of P k ,r k one obtains P = ( o° o)' r ° = PoHo a o, 

where go = S + Ro so that Bq = I 9 ° o c 2 q ° l J and therefore S C Bq = 

V 90 90 + So / 

(0,1), S C B+C^ = 1 so that P 1 = R^H'^ and n = R\H[a\. It's easy 
to prove by induction that P k = R k H' k H k and r k = R k H' k a k . Theorem 6 
implies: z T := P^r T represents the ^-minimax estimation for I 6 £(t), where 
C(t) = {Aei,ei = [1,0], A G M} if /i r ^ and £(r) = {0} otherwise. Hence 
if £ = [0, 1] then the £-minimax error is infinite. This fact reflects the non- 
causality of the model (30): since the a posteriori set of (30) is a shift of the 
convex set £f(0) = {[zq . . . z r ] : 

il(£M 2 + £il(*L t')mn 2 <,3./3>o} 

o 
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its cross section Z{r) at the instant k = r is a shift of V T (@(0)). Thus Z(t) 
is convex and unbounded implying that it recedes to infinity [61, §8] in the 
directions i ^ £(r). If h T ^ and £ = [Z,0] then the £-minimax estimation 
z T obeys 

/ w I 2 

(l, z T ) = -a T = l(z l T + (£, z T - z T f < — (31) 

since w 2 . < R~ x due to (28). Let yl, x* k denote the realization of output yk 
and state Xk derived from (28) with /& = f k , q — q* , Wk = restricted by 
(29). Let 

z l = [ z ik> z 2k]> a t be derived from (30) provided that fk = f k , 
z 2k := v k{ z ik)i Q = 1* an d w k — w t- By direct calculation z* k = x* k and 
a* k = yl so that lx* = (£, z*) with £ = [I, 0]. Hence 

(k - (/,*))»= (fr < ^ 

due to (31). Thus r h- > U,z T ) gives the online estimation of r h-> lx* with 
worst-case error ^p-- 

Minimax estimator and H 2 / filters. In [ ] a connection between set- 
membership state estimation and approach is described for linear causal 
DAEs. The authors note that the notion of informational state (X(t) in our 
notation) is shown to be intrinsic for both approaches: the mathematical 
relations between informational states of ifoo and set-membership state es- 
timation are described in [62, Lemma 6.2.]. Comparisons of set-membership 
estimators with and other widely used filters for linear DAEs are pre- 
sented in [63] provided that Fk = E. 

In [18] authors recover Kalman's recursion to LTV DAE from a deterministic 
least square fitting problem over the entire trajectory: if rank ^ k = n then 
the optimal estimation Xi\k can be found from 

Xk\k = Pk\kF'kAk~lCk-lXk-l\k-l + Pk\kH' k Rkyk, 

%o\o = Po\oH' Roy , A k l = (S^ 1 + C k Pk\kC' k ) 

Pk\l = F'kA-k-iFk + H' k R k H k , = FqSFq + H' RqH 

Corollary 5. Let r$ = F^Sq + H^R^y®. 7/rank[ j^] = n then Ik = and 

Pk Tk = Xk\k- 
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3 Numerical example 



Let us show how to use the minimax estimation in the infinite-horizon setting. 
Consider the following DAE 



o 

k 

Vi,h 

V2,k 
V3,k 
V4,k 



x l,k+l 
x 2,k+l 
x 3,k+l 

hl,k 
h-4,k 

h(i,k 



J_ 
40 
J_ 

10 



V x 3,k J V f 2 ' k J 



h-2,k 

0.005 

h 7 ,k 










x l,k 
x 2,k 
x 3,k 



+ 



91, k 

92, k 

93, k 

94, k 



(32) 
(33) 



6k 

10' "i-o 
1000, h 2 , k -- 
h,o = 0, h 7 ,k = 10fc, h 7t0 
—3 and suppose that 



where hij 
lOOfc, h 4fi 



6 

10' ,4 3,fc 

= k, h 2 ,o = 
= 0, h 8 ,k 



150k if k is odd and otherwise; = 
'- ^5,fc = Y§o ' ^5.o = 2-jjj, h 6tk = 0.05, 
,0 = 1. Also we set x^o = 1, £2,0 = 



0, /ig 



(5g, q) + ^(^Rfefi-fc, s-fe) + (S k fk, fk) < 1 



dla §{ 35(fe+l) ' TO(fe+i)}' 5 ~ dia §{^' iio}' 
= 1 so that A Q = {£ : Q%Q N £ = 0} 



where = ^-diag{^-, ^, i, S k ■ 
Note that rank = 2 and J 2 fc+i 

is nontrivial: i = (0, 0, 1) belongs to Nq. Theorem 6 implies £(&y k+l ) = 
(—00, +00) so that the a-posteriori minimax error in the direction i is in- 
finite. Thus the estimation error is unbounded in general case. Really 
(4 Qik+i r 2k+i) = so that 



|(^> Qtk+i r 2k+l — %2k+l] 



F3,2fe+l| 



Note, that any function k 1— > x 3)k satisfies (32). In this sense model (32)- 
(33) is non-causal. Since the estimation error in the direction I coincides 
with x 3j fc its evolution is unpredictable for odd k. In this sense the subspace 
{ct£, a G K 1 } in the system state space is not observable for odd k. On 
the other hand I 2k — 0. Thus R(Q 2 k) — ^ 3 and the system state space is 
observable in any direction I = R 3 due to Theorem 6. 

The dynamics of Xi,k>&i,k> \ x i,k ~ x i,k\, i — 1,3 and minimax error is 
illustrated by figures 1-2. Note, that Figure 2 demonstrates a singular case: 
for even k minimax estimation and error vanish but for odd k they gives 
nontrivial approximation. Thus one can observe some kind of oscillation of 
the estimation curve: x 32 k is near x 32k and x 32 k+i = 0. Note, that although 
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minimax error in the direction £ = (0,0,1) is infinite but (Q^k+i^^) = 0- 
Thus the corresponding minimax error curve gives an upper bound of the 
\ x 3,k ~ £3,fc| for odd k and vanishes for even k. 
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Figure 1: N = 40, state X2,k (solid), estimation x 2 ,fc (dashed); real error 
\%2,k ~ %2,k\ (dashed), minimax error (solid). 
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Figure 2: N = 40, state x^k (solid), estimation x^k (dashed); real error 
\ x 3,k — %3,k\ (dashed), minimax error (solid). 
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